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We consider a general central-field system in D dimensions and show that the division of the 
kinetic energy into radial and angular parts proceeds differently in the wavefunction picture and 
the Weyl-Wigner phase-space picture. Thus, the radial and angular kinetic energies are different 
quantities in the two pictures, containing different physical information, but the relation between 
them is well defined. We discuss this relation and illustrate its nature by examples referring to a 
free-particle and to a ground-state hydrogen atom. 
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Q \ I. INTRODUCTION 

O: 

; Phase-space representations of quantum mechanics play an increasingly important role in several branches of physics, 
' including quantum optics and atomic physics. The principal reason for this is the conceptual possibility these rep- 
resentations give for viewing the position and momentum characteristics of a quantum state in the same picture. 
^— I ' Phase space is often useful for the description of stationary states, and it has become a natural background for de- 
^ ' scribing the quantum-mechanical time evolution of wavepackets, for both matter waves and electromagnetic waves. 
Several phase-space representations have been discussed in the literature, but one of them — the so-called Weyl-Wigner 
\ representation — has come to play the role of a canonical phase-space representation, because of its simplicity . In 
accordance with this, we shall exclusively consider the Weyl-Wigner representation in the following. 

We consider this phase-space representation to be a representation in its own right. In previous work we have 
^ justified this statement by analyzing and solving the phase-space differential equations that the Wigner functions must 
^ satisfy. In particular, we have stressed that the Wigner functions may be determined directly from these equations, 
^ I without reference to wavef unctions — although it is in general easier to determine them from the wavef unctions. 
^ ^' The fact that the phase-space description is a representation in its own right makes it relevant to apply physical 
I _ intuition to the form and behavior of the Wigner functions, just as physical intuition may be applied to the form 
^ i and behavior of wavefunctions. When we do this, we discover that our understanding of quantum states becomes 
d ' enlarged, because the two types of intuition may work differently and therefore supplement each other. 
^ , In the present investigation which, for the sake of generality, is carried out in D dimensions, we consider quantum 
^ '' states referred to a center O. We focus, in particular, on the evaluation of the angular momentum and the kinetic 
^ , energy of such states. In the familiar picture based on wavefunctions, these quantities are calculated as the expectation 
values of operators. In the phase-space picture they are calculated by taking averages of dynamical phase-space 
functions with Wigner distribution functions. Performing the two calculations with care will, of course, lead to the 
. same result. Yet, a comparison between the detailed features of the two descriptions leads to some interesting and 
" " ' physically important observations. 

This was already noted in our previous work on the Wigner function for the ground state of the hydrogen atom [ 3 1 , 
in which we touched on a pedagogical dilemma which, for instance, has bothered writers of elementary textbooks 
How does one bring the fact that the angular momentum in the Bohr orbit is non-zero into accordance with the fact 
that the angular momentum in the Schrodinger picture is zero? We referred to this dilemma as the angular-momentum 
dilemma and showed that it could be resolved by noting that the mapping of the operator to phase space produces 
the phase-space function \r x — 3?i^/2 rather than just \r x pp. 

In the following, we generalize this result to D dimensions. In addition, we derive parallel but more faceted relations 
for kinetic-energy quantities, likewise in D dimensions. We discuss these results and show that the separation of kinetic 
energy into a radial and an angular part may be done in two physically meaningful ways. One is suggested by the 
form of the operators in the wave-function picture, the other by classical-like dynamical functions in the phase-space 
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picture. The relation between the two variants of radial and angular kinetic energies is tied to the Weyl correspondence 
rule and is, therefore, well defined. 

We illustrate the conceptual difference between the two types of kinetic-energy separation by two important ex- 
amples in three dimensions. One is the simplest possible time-dependent state of a free particle, the other is the 
stationary ground state of the hydrogen atom. For the first example, we find that the phase-space induced separation 
of the kinetic energy into a radial and an angular part depends on time in an intuitively simple way, whereas the 
operator-based separation is independent of time. For the ground state of the hydrogen atom, the operator-based 
separation leads to an angular kinetic energy of zero, whereas the phase-space induced separation classifies the whole 
kinetic energy as angular kinetic energy. This striking difference between the results of the two types of separation is 
well reflected in the form of the wavefunction versus the form of the Wigner function. It illustrates in a perfect way 
how the physical richness that is hidden in the simplest state of the simplest atom can only be seen by looking at the 
state from different angles. It also illustrates the intricate way in which the roots of classical mechanics are buried in 
the quantum-mechanical soil. 

The paper, which is intended to be reasonably self contained, is organized in the following way: In Sec. ^ we give 
a brief overview of hyperspherical coordinates and the central-field form of wavefunctions in D dimensions. In Sec. 
[II we define the angular momentum and perform the operator-based separation of the kinetic energy into a radial 
and an angular part. We express the result both in terms of general operators and in terms of differential operators. 
As a background for the rest of the paper we recall the salient aspects of the Weyl- Wigner transformation in Sec. IV 
. In Sec. we discuss the concept of angular momentum in the phase-space picture. We introduce the concepts of 
q {quantum) angular momentum and c (classical-like) angular momentum and discuss the relation between them. In 
Sec. we give a similar discussion of the kinetic energy in phase space and of its separation into a radial part and 
an angular part. Sees. VII and VIII are devoted to two illustrative examples in three dimensions. Sec. IX generalizes 
the examples to D dimensions. Sec. ^ is our conclusion. 



II. HYPERSPHERICAL COORDINATES 



Let r — (xi, X2, ■ ■ ■ , xd) be the position vector of a 'particle' moving in Z?-dimensional position space and 
let p — {pi,p2, ■ ■ ■ ,Pd) be its conjugate momentum. We take Xi, and hence also pi, to be Cartesian coordinates. 
In accordance with this, we introduce the hyperradius r by the relation — x\ + x\ + . . . + x^, and likewise p, 
the magnitude of the momentum, by the relation p^ = p\ 
position-space representation and write 



Pd- 



1 -^2 

Quantum mechanically, we adopt the 



P= {pi,P2,---,Pd) 



-ih 



d d d 

dxi ' dx2 ' ' dxD 



(1) 



and 



dx\ 



dxl 



dxl 



(2) 



where is the Z?-dimensional Laplacian. The kinetic energy of a quantum-mechanical particle with mass M is 
represented by the operator 



T 



P 

2M 



2M 



(3) 



The central-field Hamiltonian 



i? = — + V(r) 
2M ^ ' 



(4) 



determines the motion of the particle in a central field V{r). 

To introduce hyperspherical coordinates in position space, one writes Xi = rr]i, where the rji^s are D functions of 
D — 1 angular coordinates. Both the angles and the rji^s may be chosen in different ways, but a choice similar to the 
following one is generally used, albeit with varying notation for the angles: 



Xi — r smt/D^i sm ajj^2 ■ ■ ■ sm (^2 sinful 
X2 — r sin0£)_i shi9D-2 ■ ■ ■ sin02 cos 6*1 
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X3 = r sin sin 0d_2 ■ • • cos 62 



xd-1 = rsm&D-i cosc^i3_2 

xd — r cos(?D_i (5) 

with (0 < r < cx)), (0 < < 27r), (0 < 02 < tt), . . . , (0 < Od^i < tt). A similar representation may, of course, be set 
up for the vector p in momentum space. 

For the volume element in position space we have the expression 

dr = dxidx2 ■ ■ ■ dxo — r^^^drd^, (6) 

where the solid-angle element dCl is given by 

dn = (sin 6iD-i)^"^(sin 6*0-2)^"^ ■ . .sin 6*2 

xdeD-i...de2dei. (?) 

Integrating over all angles gives the total solid angle: 

dfl = — — - . (8) 

r(D/2) ^ ' 

Wavefunctions of the central-field problem are conveniently referred to basis functions of the form 

^p{r) = R{r)Y{Q) (9) 

where is a collective notation for the angular coordinates {9i, 62, ■ ■ ■ , do-i), and Y{p,) is a hyperspherical harmonic. 
The hyperspherical harmonics were introduced and extensively studied by Green Q and Hill [Q. They have also 
been much studied by later authors. (See, in particular, the comprehensive presentations by Sommerfcld Louck 
and Avery jl^.) The hyperspherical harmonics are eigenfunctions of the operator which represents the square 
of the total angular momentum and is defined below. 



III. ANGULAR MOMENTUM AND KINETIC ENERGY 



The angular-momentum tensor in D dimensions is defined by the operators 

Lij = XiPj - XjPi, i ^ j. (10) 
The square of the total angular momentum is 



D D 



lET:^r (11) 



2 

i=i j=i 

where the prime indicates that the double sum excludes terms for which i = j. The angular-momentum operators are 
independent of r. They merely depend upon the angular coordinates. 

We shall now separate the kinetic-energy operator T into two distinctive parts. To accomplish this in a manner 
that eases the subsequent transition to the phase-space representation, we begin by decomposing the _D-dimensional 
unit matrix 1 as follows: 

i=i ]=i 

The matrices S and T'^ are defined by the relations 

Ski = XkXi (13) 

and 
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{T^^)ki = XiSkjSij + x'^Skidu - XiXj{Ski6ij + SkjSu), 



or, 



S = 



X\ X1X2 ■ ■ ■ XiXd 
XlX\ x\ ■■■ X2XD 



XdXi XdX2 ■ ■ ■ Xj 



and, for instance. 



J.12 _ 



Adopting a dyadic notation, we may then write 





-X1X2 


• 


■ 


—X2X1 


xl 


• 


• 








• 


• 


V 





• 





Next, we note that 



and 



i=i j=i 



, S ^ r\ /r ^ 

P- -o ■P= (P- -j [- P 

r \ r / \r 



^ T^o ^ _ 1,2. 
Hence, the kinetic-energy operator may be written 

where Trad has the form 



2Mr2 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



It represents the radial kinetic energy. Tang, which represents the angular kinetic energy, is given by the operator 



(22) 



A modified expression for Tj-^d may be obtained by introducing the radial momentum Pr by the definition 

1 /V r\ 

and reahzing that 



(24) 



This gives 
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>- rad 



h\D-l){D-3) 
2M 8Mr2 ■ ^ ' 



To express Trad as a differential operator, we note from Eq. that the definition of hyperspherical coordinates 
implies that 

dr dr dxi r dxi r ^' ^"^^^ 

i—l 2—1 

Using this, and the commutation relations between the components of r and p, turns the expression (^l]) into the 
following forms 

™''~ 2MrD-idr dr 
f D-l d 



2M \dr^ r dr 

These are familiar expressions. 



(27) 



IV. THE WEYL-WIGNER TRANSFORMATION 

Let ip{r) be a normalized position-space wavefunction, 

'iP*{r)i;{r)dr = 1. (28) 

Further, let A be some operator acting on ■0. The expectation value of A in the state '0 is then 

(i) = / ij*{r)Aij{r)dr. (29) 



When A is Hermitian, (A) is real. 

Another way to evaluate the expectation value (|9|) is by introducing the 2D-dimensional (r,p) phase space, and 
then use the expression 

(i) = / I a[r,p)W{r,p)drdp. (30) 



Here, W{r^p) is the Wigner function corresponding to the wavefunction tp{r)^ and a{r,p) is the dynamical phase- 
space function corresponding to the operator A. For a Hermitian operator a{r,p) is real. The Wigner function 
W{r,p) is always real, but may take negative values. Its integral over phase space is, however, always equal to 1, 



W{r,p)drdp^l. (31) 
It has the particle density p{r) and the momentum density H(p) as marginal densities: 

p(r) = j W{r,p)dp, (32a) 
n(p) = j W{r,p)dr. (32b) 
The Wigner function is defined as follows @[6|-[l|0,| 



D 



nh J 

^e-'Pr'/f'dr'. (33) 



mr-,p)=(^j I nr-r'/2)^{r + r'/2) 
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If the operator A is of the form F(f) + G(p), then a(r,p) will be simply F(r) + G(p). Otherwise, the noncommutativity 
between r and p will come into play. The transformation involved is the Weyl transformation |]l5|-[l8|,|l|,^ . It may, for 
instance, be represented in the following form 



a{r,p) = / (r + r' /2\A\r - r' /2)e-'P-''' ^'''dr' , 



(34) 



in which the matrix element of A is defined with respect to two eigenstates of the position-vector operator, corre- 
sponding to the eigenvalues r + r' /2 and r — r'/2, respectively. 

Let C = AB. The dynamical phase-space function corresponding to the operator C is then given by the star 
product 



c(r,p) = a(r,p) * b{r,p), 



(35) 



where 



c(r,p) = exp 



ifi f d d d d 
2 \dri dp2 dpi dr2 
a{r,p)b{r,p). 



(36) 



Here, the subscript 1 on a differential operator indicates that this operator acts only on the first function in the 
product a{r , p)b{r , p) . Similarly, the subscript 2 is used with operators that only act on the second function in the 
product. 



V. ANGULAR MOMENTUM IN THE PHASE-SPACE PICTURE 

We denote the Weyl transform of the operator Lij by Aij. It is simply given by 

Aij = XiPj - XjPt. (37) 
In analogy with Eq. ( pi] ) we introduce the dynamical phase-space function 

^' = ^EE'4- (38) 
i=l J=l 

This function is, however, not the Weyl transform of L^. 

To determine the actual Weyl transform of , we first determine the Weyl transform of the operator Lfj . We 
do this by invoking the relation (^), with a and b both equal to Ay . It is found that only the three first terms in 
the expansion of the exponential operator contribute to the result. The Weyl transform of Lfj is thus found to be 
A^j — \'h^- We express the result as the mapping 

1% ^ Kl - \n\ (39) 



2 

From this, we get for the square of the total angular momentum: 



L^^A^-^^n\ (40) 

This is the generalisation of the result for D = 3 that we gave in the Introduction. 

Hence we have, for a state described by the wavefunction tp{r) and the Wigner function W{r,p), that 



K^W{r,p)drdp- ^^^~^\ ''. (41) 
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To put this and the other relations above into perspective, let us distinguish between two equally well defined angular 
momenta. One of these goes naturally with the wavefunction description. The other goes naturally with the phase- 
space description. The first angular momentum, which we shall call the q {quantum) angular momentum, is defined 
by the operators Lij and . The other angular momentum we shall call the c (classical-like) angular momentum. It 
is defined by the dynamical phase-space functions and A^. The two angular momenta are connected by relations 
like that of Eq. (^), but the kinds of intuition one may attach to them are quite different. 

Thus, the q angular momentum is primarily an algebraic, or group-theoretical concept. The Lij operators are 
generators of infinitesimal rotations, and their eigenvalues describe the possible behavior of a given state under 
rotations. An s-state in a three-dimensional world, with its zero-eigenvalue of L^, is for instance invariant under a 
rotation about any axis through the origin. This is what the eigenvalue of tells us. The position-space wavefunction 
'ijj{r) is independent of angles, and hence it predicts all directions of r to be equally probable. If we prefer to describe 
the state by its momentum-space wavefunction 

m = (2^) ' I tP{r)e-'P-^/''dr, (42) 

then this wavefunction is independent of angles about the origin of momentum space. Hence, all directions of p are 
also equally probable. But one cannot talk about a coupling between the directions of r and p in this description. 

The expression for the c angular momentum is, on the other hand, just the classical angular-momentum expression 
for particle motion relative to the center O. Thus, it provides a measure of the relative direction of r and p with 
respect to O. This directional correlation, which is not directly apparent from the wavefunction, is manifestly present 
in the Wigner function. It leads to a non-zero c angular momentum even if the q angular momentum vanishes. This 
is in complete accordance with the general relation (^) . 



VI. KINETIC ENERGY IN THE PHASE-SPACE PICTURE 

The Weyl transform of the kinetic-energy operator T given by the expression (H), is simply 

T"^^. (43) 
To separate it into a classical-like radial part, T^ad' ^ classical-like angular part, T^„g, we note that the resolution 



of the identity matrix given by Eq. (12) also holds in phase space. Hence, we may also write 



2EE— 



But now the components oip commute with the components of r, and the analogue of Eq. (po| ) becomes 

rpC rpC • rjiC 

rad ' ang 

1 /p-r\2 A2 



2M \ r J 2Mr2 

with A^ given by Eq. (^). The phase-space version of the radial kinetic energy is consequently 

1 (p-r^ 
2M V r 
(pcosu)^ 
~ 2M 



(45) 



(46) 



The last expression is obtained by putting r p = rpcosu, where u is the angle between r and p. The angular kinetic 
energy is 

A2 

= (P^^^")' (47) 
2M ■ ^ ' 
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The validity of the last expression is simplest verified by noting that — (pcosw)^ = (psinu)^. 

It is important to note that T^^^ is not the Weyl transform of Trad, nor is T^„g the Weyl transform of Tang- To 

determine the actual Weyl transforms of Trad and Tang requires repeated applications of the relation ( ^6| ) to the 
expressions ( ^l| ) and (]2^). We find, after some tedious algebra: 

~ {D-l){D-2)h' 

Irad' >-L rad + (48) 



and 



f (49) 

J-ang' > J- ang SMr^ " ' 

The relation (43) between the radial kinetic energies might also have been obtained from the relation (]25|) by 



exploiting the fo 



and 



lowing interesting relations: 



^ (50) 



^2 /r-p\2 {D-l)h'^ 



r 



Pt^i — ) 4,2 ^ (51) 



where pr- is the operator defined by Eq. (^). These relations may again be derived by repeated application of the 
expression (^) for the star product. 

We shall now apply the results of this and the previous section to two important examples. 

VII. A FREE-PARTICLE STATE 

In this section wc consider a minimum- uncertainty state in three dimensions. The wavcfunction for such a state is 

^{r,0)^^—j e-^"'-, (52) 
at a chosen initial time t = 0. The corresponding Wigner function has the form 

W{r,p,0)^^e-'^'^'-p'^-'^\ (53) 

Let us write down the Wigner function at a later time t under the assumption that we are dealing with a free 
particle. We then know that the Wigner function will evolve in time in exactly the same way as a classical phase- 
space distribution, that is, it will develop according to the classical Liouville equation Thus, the expression 
for W(r,p,t) may be obtained from that for W{r,p,0) by simply replacing r by r — tp/M. In this way we get 

W{r,p, t) = -L^^-o?r--({tlrf+X)v-i^-n-+2i,tiry.pin^ (54) 

where 

Since r p — rpcosu, where u is the angle between the vectors r and p, the Wigner function (|5j) only depends on the 
three coordinates r, p and u. It is independent of the three Euler angles that determine the orientation of the {r,p) 
cross in phase space. This independence of the Euler angles expresses the overall rotational invariance of the state, 
that is, it signifies that the oangular momentum is zero at all times. However, evaluating the average of with the 
Wigner functions ( p3| ) and ([sj) gives as it should according to the relation (p|). 

We note that all angles between r and p are equally probable at t = 0, whereas small values of u are favored for 
large values of t. Thus, the correlation between r and p depends on t. Such a statement would be impossible to 
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defend in the wavefunction picture. However, in the phase-space picture it is very meaningful. For at i = 0, the 
Wigner function (|5j) is concentrated about the origin of phase space. But in the course of time, phase space points 
with non-zero values of p will move to points with larger values of r, with r and p becoming more and more parallel. 
But this just amounts to small values of u being favored for large values of t, as in the expression js^). 

This behavior is also reflected in the expectation values of the radial and angular parts of the kinetic energy as a 
function of time, but only in the phase-space picture. It is readily verified that the total kinetic energy associated 
with the wavefunction ( p2| ) has the value 

It is, of course, independent of t. And since ipir, t) represents an s state, the expectation value of the fang of Eq. ( [2^ ) 
is zero at all times. Thus we have: 

IT 

^ ™''' ~ 4 M ' 

{fang) = 0, (57) 

for all t. 

In the phase-space picture, we must average the dynamical phase-space functions T'^, T^^^ and T^^g with the 
phase-space function(^). We denote the resulting energies by T, Trad and 7^„g, respectively: 

1 t t .2, 



and 



We have, of course, that 



'^^2mJ J P^i'-^P^drdp, (58) 
%ad = ^ J J{pcosufW{r,p)drdp, (59) 

%ng J ipsinu)^Wir,p)drdp, (60) 



T={T). (61) 

Rather than evaluating the expressions (^9|) and (|6^) by direct integration, we may proceed in the following 
analytical way: According to Eqs. (^) and (^9|), with D = 3, T^^id ^^^'^ '^ang are the Weyl transforms of the operators 

^2 

Trad = Trad ~ (62) 



and 



Trad - Tang + ^j^^2 ' (^^) 



respectively. Hence, we merely have to modify the values in (57) with the expectation value of fi^ /AMr^ to get the 
values of Trad and Tang- Evaluating this value by averaging h /AMr'^ with the Wigner function (Q) gives 

/^U^ ' . (64) 
\4Afr2/ 2M (t/T)2-M ^ ' 



Hence, we get 



T - " 

■'■rad — 



M V4 2((t/T)2 + l) 
_ a2^2 ^ 

M 2((t/T)2 + l)- 

Thus, the classical-like kinetic energy is purely radial for very large values of t. However, at t = only one third of 
the kinetic energy is radial, corresponding to the value (l/4)(a2?i^/M). The angular kinetic energy is twice as large. 
This reflects the fact that r and p are entirely uncorrelated at t = 0, so that p is twice as likely to be perpendicular 
to r as being parallel to r. 

The energies T, Trad and Tang are shown graphically in Fig. |l|, as functions of t. 
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VIII. THE HYDROGEN ATOM 



Having studied a free-particle case, we shall next consider a bound-state case, namely, the ground state of the 
three-dimensional hydrogen atom. The wavcfunction is now 



^{r) ^ J ^e-^/'^\ (66) 



Trap 



where ao is the Bohr radius, 







h"^ 47ren 



M is the electron mass (we treat the nucleus as being infinitely heavy), and e is the magnitude of the elementary 
charge. The state considered is a stationary state, and the Wigner function is accordingly independent of time. As in 
the previous example, we are dealing with an s state, so the Wigner function is again independent of the Euler angles 
that determine the orientation of the (r,p) cross, but it does depend on the angle u between r and p, and in fact in a 
more complicated manner than in the expression (|5^). In addition, the Wigner function for the hydrogen atom takes 
both positive and negative values, whereas the Wigner function (|5j) is non-negative at all times. 

We have made a detailed study of the Wigner function for the hydrogen atom in Our results were, inter alia, 
presented as a series of contour maps for different u values. These maps showed that the Wigner function is everywhere 
positive in what we called the dominant subspace, that is, the part of phase space in which r and p are perpendicular 
{u = vr/2). It is, in particular, large in the part of this subspace obtained by putting r — gq and p — h/ap. This is the 
region of phase space to which the ground-state motion was restricted in early quantum mechanics Il9|, since a Bohr 
orbit (in position space) is just a circle with radius aq, in which the electron is supposed to move with the constant 
momentum h/ap. For other angles than u = n/2, the Wigner function develops negative regions. It is, in particular, 
strongly oscillating around the value zero for small angles and angles approaching tt. 

We shall now see that this pronounced correlation between the directions of r and p is strongly reflected in the 
partitioning of the kinetic energy. The kinetic energy associated with the wavcfunction ( p6[ ) has the value 

and since ip{r) represents an s state we have, in analogy with the expressions 



(Trad) — 



2Mal ' 



{fang) - 0. (69) 

The relations ( |6^ and (^) still hold, and it is readily found that 



2Mal 



(70) 



The analogue of (|65| ) becomes therefore 



%ad — 0, 



^""^ " 2Mal ■ ^'^^^ 

This is a remarkable result. For it shows that, in the phase-space representation, the radial kinetic energy vanishes. 
The kinetic energy is purely angular. This is of course in complete harmony with a picture in which the electron 
primarily revolves around the nucleus rather than moving in the radial direction. 

The cases studied in this section and the previous one have both been for s states. The expressions we have derived 
in the first six sections are, however, valid for any state independent of its angular momentum, but the interesting 
effects are most pronounced for the s states. With a minor exception, the derived expressions are also valid for any 
dimension D. A few comments concerning dimensions different from three are, therefore, in order. 
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IX. ARBITRARY DIMENSIONS 



The minor exception mentioned above has to do with zero- angular- momentum states for D = 2. The relations ( |4^ ) 
and (^9|) suggest that the separation of the kinetic energy into a radial part and an angular part is invariant under 
the Weyl transformation for D = 2. This is, however, not quite true. One must be aware that taking expectation 
values with the expressions ( ^8| ) and (^) involves taking expectation values of l/r^ , and such expectation values are 
undefined for wavefunctions that stay finite at r = in a two-dimensional world. This is because the volume element 
only contains r to the first power for D = 2. Hence, one cannot exploit operator relations like those of Eqs. ( |6^ ) 
and (16^ ) for zero-angular-momentum states. This, however, does not reduce the significance of the dynamical phase- 
space functions T^^^ and T^ng- It just implies that the radial and angular kinetic energies Trad and Tang associated 
with them must be evaluated directly by averaging with the Wigner function, using the two-dimensional version of 
the expressions (^) and (|^. For the two-dimensional equivalent of the Wigner function (|5j) this produces a t 
dependence of the kinetic energies similar to that in Fig. with the difference that the kinetic energy is equally 
distributed on its radial and angular components at t = 0. This difference was to be expected since there is only one 
perpendicular direction to r in two dimensions. 

Concerning the hydrogen atom, it is interesting to note that the above conclusions for the ground state of the 
hydrogen atom in three dimensions remain valid for other values of D. Thus, the kinetic energy is purely radial in 
the wavefunction picture, but purely angular in the phase-space picture. The ground-state wavefunction for the D 
dimensional hydrogen atom has the general form 



*'-'^\/i »y(n„/2)'...,..(2„„)! -'"^'-° 



where Sjj is the total solid angle (ph and 



The ground-state kinetic energy is 



D -I 

no = ■ (73) 



and for D > 3 exactly the same value is found for the expectation value of the D dependent term in (48) and (|4£ 
This confirms that the classical-like kinetic energy is, in fact, purely angular for D > 3. 

For D = 2,we can again not draw on expressions like those of Eqs. ( p2[ ) and (|63[). We must perform the phase-space 
integrations directly. Doing so shows that also for D = 2, the classical-like kinetic energy is purely angular. 

The said integrations over phase space are far from simple to perform, because the Wigner function for the hydrogen 
atom cannot be evaluated analytically A practical procedure is to expand the Is wavefunction on a set of 

Gaussians. For a linear combination of N Gaussians, the Wigner function may be determined analytically. The 
values of Tad and Tng may then be calculated by combining analytical and numerical integrations. For = 1, 
we have a Wigner function similar to that of Eq. (^3|) for a free particle, but in two dimensions only, leading to 
Tad = Tang- As N incrcascs, the contribution from Tad is found to decrease, converging to zero for large values of 
N. This is shown graphically in Fig. |2[ 

The information in Fig. ^ is not merely numerical. It serves as yet another demonstration of the physical uniqueness 
of the Coulomb potential. For instead of considering the used wavefunctions to be approximate solutions for the 
Coulomb potential, we may consider them to be exact solutions for a different potential. As A'^ increases, this potential 
becomes more and more Coulomb like. This suggests that only for the Coulomb potential is the classical-like kinetic 
energy purely angular. 

The linear combinations of Gaussians used to prepare Fig. |2| were determined by the variational principle, following 
the prescription given in reference pO[ |. In that work, we made explicit studies of the Wigner function for the D 
dimensional hydrogen atom and presented contour curves for selected values of D. These contour maps show, inter 
alia, that the oscillations between negative and positive values of the Wigner function become weaker and weaker for 
higher D values. We may therefore say that the phase-space distributions become more classical as the dimensionality 
D increase. 
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X. DISCUSSION 



The problem of separating the kinetic energy of a particle moving in a central field into a radial part and an angular 

part has not attracted much attention in the past. It has been tacitly assumed that the separation offered by the 
wavefunction picture was the only sensible one. In the present paper we have challenged this view by putting focus 
upon the separation offered by the phase-space picture. This separation is equally well defined. In contrast to the 
former, it throws much light on the correlation between the directions of the position vector r and the momentum 
vector p. There is no way of discussing this correlation in the wavefunction picture. 

The most amazing result of our analysis is probably our finding that the kinetic energy in the ground state of 
the hydrogen atom is purely radial in the wavefunction picture, and purely angular in the phase-space picture. This 
finding, which relates to arbitrary dimensions, reflects the great difference in the kind of intuition one may apply 
within the two pictures. This difference is also well illustrated by focusing upon the angular momentum. In the 
wavefunction picture, the angular momentum primarily refers to the behavior of the state under rotations. In the 
phase-space picture it also refers to the correlation between the directions of r and p. 

To properly understand a quantum system one must look at it from different angles. The present paper shows once 
again that the phase-space picture is a fruitful supplement to the wavefunction picture. 
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FIG. 1. Time dependence of the classical-like kinetic energy of a free particle described by the Wigner function (JSJ). t 
measured in units of r = M/a^h, energies in units of e = a^h? /M. The kinetic energy T is separated into its radial part 7^ 
and its angular part Tang- 
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FIG. 2. With the ground-state wavefunction of the two-dimensional hydrogen atom approximated by a hnear combination 
of N Gaussians, the figure shows the A^-dependence of the calculated classical-like kinetic energy T, as well as its radial and 
angular parts, Trad and Ta„g. Energies are measured in units of eo = /MaQ. 
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